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In this paper, we calculate the effective gravitational mass of Ayon-Beato and Garcia regular 
(non-singular) static spherically symmetric asymptotically Minkowskian metric that is a solution 
to Einstein’s equations coupled with a nonlinear electromagnetic field. The effective gravitational 
mass is negative, zero, or positive that depends on the ratio of magnitude of electric charge to the 
ADM mass and the ratio of the radial distance to the ADM mass. As expected, at large value of 
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I. INTRODUCTION 


In the year 1984, Cohen and de Felice [T] defined effective gravitational mass of a metric at 


any point in a given space-time at a radial distance r as the energy contained inside the region 
of r = constant surface, and used this to explain repulsive nature of t im elike singularities. 
They used Komar energy formula and obtained effective gravitational mass of the Kerr- 
Newman metric (characterized by mass M, electric charge q, and rotation parameter a): 



( 1 ) 


For large values of rotation and/or charge parameters (that is, when these parameters dom¬ 
inate over the mass parameter), the effective gravitational mass becomes negative. For the 
Reissner-Nordstrom metric (a = 0) in the above result, the effective gravitation mass is 


H(r) = M 


( 2 ) 


2 

Obviously, /i(r) < 0 for r < |^ and this explains repulsive gravitational effect on an electri¬ 


cally neutral test particle in the gravitational held described by the Reissner-Nordstrom met¬ 
ric. This fascinating interpretation of energy content in a space-time attracted researchers; 
however, they soon realized that there is no unique adequate formalism for energy calculation 
and this had been a difficult problem since the appearance of general theory of relativity. 
For example, Komar definition cannot be applied to non-static space-times. 

Though Einstein’s theory of gravity is a well governing theory of space, time, and grav¬ 
itation and it is very well experimentally verified, and is, so far, proven to be the most 
successful theory to describe gravity at small as well as large scales, it does possess some is¬ 
sues. The energy-momentum localization problem is one such important issue in the context 
of General Relativity. Energy, momentum, and angular momentum are important conserved 
quantities in Minkowski space-time (when gravitational held is negligibly small.) These have 
significant role as they provide the hrst integrals of motion to solve otherwise unmanage¬ 
able physical problems. However, till today, we do not have a general definition to describe 
energy-momentum distribution in curved space-times (i.e., in presence of gravitational held) 
and several difficulties arouse in this direction. 

All the attempts to hnding this resulted in a large number of different energy-momentum 
complexes. To solve this problem, Einstein formulated the energy-momentum local conser¬ 
vation law (see in 0 )- Then, many physicists including, Landau-Lifshitz[3], Weinberg^, 
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Papapetrou[5j, and Bergman and Thomson [6 J proposed various definitions of energy- 
momentum distribution. But, these definitions of energy-momentum complexes were coordi¬ 
nate dependent, that is, they give meaningful results only when the calculations were done in 
quasi-Cartesian coordinates. It was suspected that a plethora of different energy-momentum 
complexes would give acceptable total energy and momentum for isolated systems (i.e., 
asymptotically flat space-times); however, they would produce different and hence meaning¬ 
less energy-momentum distributions even in asymptotically flat (Minkowskian) spacetime 
and would not give any meaningful result at all for asymptotically non-Minkowskian space¬ 
times. Virbhadra’s and his collaborators j3, |S] seminal work shook this prevailing prejudice 
and they explicitly showed that different energy-momentum complexes give the same results 
for many spacetimes (asymptotically Minkowskian as well as non-Minkowskian). Later, 
Virbhadra and Rosen[9] (the most famous collaborator of Albert Einstein) studied Einstein- 
Rosen gravitational waves and for this metric as well different complexes gave same and 
reasonable results. Further, attracted by these encouraging results, many researchers [TUI [T2] 
studied several other metrics and obtained very useful results. 

In General Theory of Relativity, majority of well-known exact black holes solutions to the 
Einstein’s field equations came up with the existence of space-time (curvature) singularities 
where general relativity theory breaks down. Therefore, it is very desirable to have regular 
(nonsingular) solutions of Einstein’s equations. To avoid this black hole curvature singularity 
aspect, researchers came up with a large number of regular black hole models which were 
referred to as “Bardeen black holes” [T3] because Bardeen was the first to obtain a regular 
black hole solution. But, the problem with all of them was that none of them were an exact 
solution to the Einstein’s Field Equations. Moreover, there are no known physical sources 
pertaining to any of those regular black hole solutions. 

Later in 1998, Ayon-Beato and Garcia found out the first singularity free exact black hole 
solution[Tl] to the Einstein’s field equations whose source is a non-linear electrodynamic field 
coupled to gravity and the metric does not possess any curvature singularity. In this paper, 
we calculate effective gravitational mass of this metric and analyze the results. We use 
Mpllcr’s [2j energy-momentum expression which permits to calculate the energy-momentum 
distribution in an arbitrary coordinate system. As usual in general relativity papers, we too 
use the geometrized units (G = 1, c = 1). We use Mathematica [H] software for plots. 
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II. THE AYON-BEATO AND GARCIA SOLUTION 


The Ayon-Beato and Garcia (hereafter referred to as AG) static spherically symmetric 
asymptotically Minkowskian solution to Einstein’s equations with a source of nonlinear 
electrodynamic field is given by the line-element 


ds 2 = fddt 2 — adr 2 — r 2 {dd 2 + sin 2 6d(ft 2 ) 


(3) 


and the electric field 


where 


4 i r 2 — 5 q 2 15 M 

E = qr 4 ( —- q — + 


(r 2 + g 2 ) 4 2 ( r 2 _)_ g 2 ) 7 / 2 J 


P=- = l 


2 Mr 2 


+ 


q 2 r 2 


(4) 


(5) 


a (r 2 T O' 2 )i (r 2 + g 2 ) 2 ’ 

Ayon-Beato and Garcia na showed that their solution corresponds to charged non-singular 
black holes for \q\ < approximately 0.6 M and otherwise the solutions are though regular, 
there are no horizons. 

The asymptotic behavior of the AG solution is given by 

2 M 4 M 2 -q 2 „( 1 

a = 1 + -+- jt-E + 0 [ — 

r 

2 M , q 2 
r 2 


f3 = 1 -— + ^+ 0 (^ 1 ; an, 4 


r 

E = 


( 6 ) 


Thus the A-G solution behaves asymptotically as the Reissner-Nordstrom solution to 
Einstein-Maxwell equations. The symbol M represents the ADM mass and q stands for 
the electric charge. The behavior of the above functions near r = 0: 


a 

V 

E 


1 + r 2 
1 — r 2 

/ 15|g| 
V 2« 6 


- y) r “ + 0 ■ 


(7) 


Thus, the metric behaves as de Sitter near r = 0. As we will need the metric in quasi- 
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Cartesian coordinate system as well, applying the coordinate transformation 

x = r sin 6 cos 0 , 
y = r sin 9 sin 0, 

z = rcosO , ( 8 ) 


the line element gets the form, 

ds 2 = (3dt 2 — ('yx 2 + l)dx 2 — ('-yy 2 + l)dy 2 — ('yz 2 + l)dz 2 — 2^(xydxdy+yzdydz + zxdzdx), (9) 


where 


A-l 


7 = 


x 2 + y 2 + z 1 


( 10 ) 


III. M0LLER ENERGY-MOMENTUM PRESCRIPTION 


The Mpllcr energy-momentum complex is given by [2j 


% kl — _ Y kl l 

which satisfies the local conservation laws as 


dx k 

where the antisymmetric Mpllcr super-potential xf is 


( 11 ) 


( 12 ) 


x? = -x? = v^5 [9in, m - g km g"‘ ■ (13) 

The energy and momentum components are now given by 


Pi 



dx x dx 2 dx 3 . 


(14) 


Pi stands for momentum components Pi, P 2 , P 3 , and Pq is the energy. Following Cohen 
and de Felice [lj, we interpret the energy function as the effective gravitational mass /i(r) 
of the metric. We apply Gauss’ theorem in above equation and then the total energy and 
momentum components take the new form: 


Pi = ~ 
87T 



Xi S n ,5 dS , (15) 

where ns is the outward unit normal vector over an infinitesimal surface element dS. 
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IV. CALCULATIONS 


M0llcr claimed that his energy-momentum complex can be used in any coordinate system 
and would produce the same results. In order to obtain the effective gravitational mass of 
the AG metric and to verify Mpller’s claim, we first perform calculations in Schwarzschild 
coordinates {t,r, 0,(f)} and then in quasi-Cartesian coordinates {t,x,y,z}. 


A. Calculations in Schwarzschild coordinates 


In Schwarzschild coordinates, the determinant of the covariant metric tensor for the 
AG metric given by Eq. (|3]) is 

g = — r 4 sin 2 0. (16) 

As the metric is represented by a diagonal matrix, g lk = A-V values of indices i and k. The 
only non-vanishing component of xT which is needed for the energy calculation is 


Xo 1 = —Xo U = 2r a sin0 


,10 


M 


2 q 2 


(r 2 + g 2 )l * (r 2 + g 2 ) 3 


(17) 


while all the other components of \ kl vanishes. Using the above expression in equation (15), 
the energy distribution is given by 


Ss = r 3 


M- 


2 q 2 


i 


- q 


(18) 


(r 2 + <j 2 )5 ‘ (i- 2 + ? 2 ) 3 . 

The subscript S emphasizes that calculations have been performed using Schwarzschild 

coordinates {t, r, Thus, we find the same result for Mpllcr energy in Schwarzschild 

coordinates as obtained by Yang et al. in [16J. (The same research team also obtained 
Mpller energy in Schwarzschild coordinates for a more general metric m-) Now we calculate 
3 momentum components and all are found to be zero, i.e., 


P r = Pe = P<f> = 0, 

and this is desired result as the AG metric is static. 


(19) 


B. Calculations in quasi-Cartesian Coordinates 

The desired covariant components of fundamental metric tensor g& is obtained from the 
equation ([9]): 
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o \ 


9ik 


( p 0 

0 —7X 2 — 1 
0 -xy'y -7 y 

\ 0 — xz'y -yz-f — 7 z l — 1 ) 

The determinant of this metric tensor is, 


0 

-27/7 
■ 2 1 


-xz^f 

-yz 7 

v 2 


( 20 ) 


9 \9ik\ 1 - 


( 21 ) 


The components of g lk are given by 




9' k = 


0 0 0 

0 —/3 (7y 2 + 7Z 2 + 1 ) xy/3 7 xz/3 7 

0 —(3 (7X 2 + 7Z 2 + 1 ) yzPj 


\ 


( 22 ) 


Y 0 xzft 7 yzPx —/3 (7X 2 + 7y 2 + 1 ) / 

The components of the super-potential xf in quasi-Cartesian coordinates are obtained here: 


01 

dp 


Xo 

dx ’ 


02 

dp 


Xo 

dy ’ 


03 

dp 


Xo 

~ dz ’ 


, 01 

„ 02 

„ 03 

Xi 

= Xi = 

Xl 

01 

02 

03 

X 2 

= X 2 = 

X 2 

„ 01 

„ 02 

, 03 

X 3 

~ X 3 — 

x 3 


(23) 


Using these x{ k info the equation (15), we obtain the energy-momentum components in 
quasi-Cartesian coordinates, 


= r 3 


M 


2 q 2 
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q 


(r 2 + g 2 ) I (r 2 + q 2 ) 3 


(24) 


The subscript QC emphasizes that calculations have been performed using quasi- 
Cartesian coordinates {t,x,y,z}. All the three momentum components are found to be 
zero, i.e., 

P x = P y = P z = 0, (25) 
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M 0 llcr energy in Schwarzschild coordinates obtained by Yang et al [16] and our results, 


given by equations (18), (19), (24), and (25) for energy and momentum distributions in 
AG space-time, support the claim of Mpllcr that his definition of energy and momentum 
is applicable to both Schwarzschild coordinates {£, r, 9, cp} and quasi-Cartesian coordinates 
{t,x,y,z}. 


C. Effective gravitational mass of the AG metric 


In the previous sub-section, we showed that Mqllcr’s definition of energy-momentum pro¬ 
duces coordinate-independent results. Now, following Cohen and de Fclice[T], we interpret 
the energy distribution as the effective gravitational mass /i(r), and then thoroughly analyze 
the results and explain its importance. 

The effective gravitational mass of the AG metric is 


ji(r) = r ' 3 


M- 


2 q 2 


q 


q 


(r 2 + q 2 )% ' (r 2 + q 2 ) 3 

The asymptotic behavior of the effective gravitational mass is 

9 Mq 2 4g 4 75 Mq A 9 q e „ ( 1 


H(r) — M — — 
r 


+ + 


A + O v 


(26) 


(27) 


2^2 rjr> 3 iy»0 

Thus the asymptotic value of the effective gravitational mass is the ADM mass M. It is 
also clear that at large r, the effective gravitational mass of the AG metric becomes the 
effective gravitational mass of the Reissner-Nordstrom metric as obtained by Cohen and de 
Felice [IJ using Komar energy and by Virbhadra[7] using Mqllcr’s complex, q — 0 in the 
above equation gives the effective gravitational mass of the Schwarzschild metric. 

The behavior of the effective gravitational mass near r = 0 is 


5 (6M\q | 

W q A ) v q 6 



+ O (r 7 ) . 


(28) 


Thus, depending on the values of ^ and jj, the effective gravitational mass can be neg¬ 
ative, zero, or positive. We now make numerous plots to see the behavior of the effective 
gravitational mass at large as well as small radial distances and also to see the facinating 
role of jj on the effective gravitational mass. 

In Figure 1, we plot the ratio of the effective gravitational mass to the ADM mass 
against the ratio of the radial distance to the ADM mass for different values of \q\/M. The 














Black Hole No horizon 




FIG. 1: (color online). In the figure on left side, the ratio of the effective gravitational mass /x to 
the ADM mass M is plotted against the ratio of the radial distance r to the ADM mass M for 
Q/M = 0.1 (blue), 0.2(green), 0.3(orange), 0.4(red), and 0.5(black). Further, the figure on right side, 
the ratio of the effective gravitational mass /x to the ADM mass M is plotted against the ratio of the 
radial distance r to the ADM mass M for Q/M = l(blue), 2(green), 3(orange), 4(red), and5(black). 
Plots on left and right sides are, respectively, for regular charged massive black holes and for regular 
charges massive objects with no horizons. 


No horizon No horizon 




FIG. 2: (color online). The ratio of the effective gravitational mass /x to the ADM mass 
M is plotted against the ratio of the radial distance r to the ADM mass M for \q\/M = 
1 (blue), 2(green), 3(orange), 4(red), and 5(black) in the vicinity of the center. 

asymptotic value for each case is 1 showing that the total mass of the metric is the ADM 
mass M. For any fixed value of the effective gravitational mass is smaller for higher value 
of jj. Thus, the electric charge contributes negatively to the effect gravitational mass. As, 
for black holes, jj is lower than for charged regular objects with no horizon, for any fixed 
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FIG. 3: (color online). The ratio of the effective gravitational mass to the ADM mass, i.e., fi/M 
is plotted against the ratio of the radial distance to the ADM mass, i.e., r/M and the ratio of the 
charge to the ADM mass \q\/M. The figure on left side is a surface plot whereas the one on right 
side is the corresponding density plot. This shows behavior of p/M at large distances. 

value of jj (but outside the even horizon), black holes have higher effective gravitational 
masses compared to those without event horizon. In Figure 3, the same effects are shown 
more clearly through surface and density plots. 

In Figure 2, We plot the same quantities as we plotted in Figure 2; however, we plot for 
no horizon cases only and near the center (r = 0). Near the center, unlike the cases of large 
radial distances, higher value of [i/M for lower value | q\/M does not always hold. In fact, 
near r = 0, p/M could be negative, zero, or positive and that is determined by both \q\/M 
and r/M. The same quantities are plotted in figures 4 and 5 as surface and density plots 
which show with more clarity the variation of the effective mass to the ADM mass ratio as 
a function of the radial distance to the ADM mass as well as the ratio of the absolute value 
of electric charge to the ADM mass. 
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No horizon 



Surface plot 



FIG. 4: (color online). The ratio of the effective gravitational mass to the ADM mass, i.e., q/M 
is plotted against the ratio of the radial distance to the ADM mass, i.e., r/M, and the ratio of the 
charge to the ADM mass \q\/M. The figure on left side is a surface plot whereas the one on right 
side is the corresponding density plot. This plot shows the behavior of /i /M close to the center. 

V. SUMMARY 

The AG solution is a static spherically symmetric and asymptotically Minkowskian reg¬ 
ular solution of Einstein’s equations coupled with a nonlinear electrodynamic field. The 
solution behaves as Reissner-Nordstrom solution at a large distance and de sitter as r ap¬ 
proaches zero. For |g| < approximately 0.6 M, this represents black holes and for larger 
values of \q\/m there is no event horizon. 

We calculated the effective gravitational mass of the AG metric and extensively ana¬ 
lyzed the results. At large distances, the role of the electric charge \q\ is to decrease the 
effective gravitational mass. However for small value of radial distances, near r = 0, the 
effective gravitational mass dependence on the electric charge is not necessarily similar and 
the effective gravitational mass could be positive, zero, or negative depending on the radial 
distance to the ADM mass ratio and relative magnitude of absolute value of electric charge 
compared to the ADM mass. As we already discussed, Yang et al H7J obtained Mqller 
energy in Schwarzschild coordinates. Our investigation in this work evidences that Mpller’s 
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No horizon 



FIG. 5: (color online). The ratio of the effective gravitational mass to the ADM mass, i.e., [i/M 
is plotted against the ratio of the radial distance to the ADM mass, i.e., r/M and the ratio of 
the charge to the ADM mass \q\/M. The figure on left side is a surface plot whereas the one on 
right side is the corresponding density plot. This plot exhibits the behavior of /j/Af very close the 
center. 

prescription to find energy and momentum is unaffected by the choice of different coordinate 
systems. 

The study of effective gravitational mass that an electrically neutral test particle experi¬ 
ences predict many physical effects even before calculations. A metric exhibiting negative ef¬ 
fective gravitational mass has repulsive effects not only to timelike but also to null geodesics, 
and hence also affect gravitational lensing phenomenon. Virbhadra’s pioneer research [T8] in 
gravitational lensing as a tool to propose astronomical test to the unproven cosmic censorship 
was influenced by his works and analysis of the energy distribution (effective gravitational 
mass) in spacetimes. The knowledge of the effective gravitational mass of AB metric would 
give more physical insights about the spacetime which might have significant applications 
to relativistic astrophysics. 
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